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Abstract. Gauge invariance in discrete dynamical systems and its con- 
nection with quantization are considered. For a complete description of 
gauge symmetries of a system we construct explicitly a class of groups 
unifying in a natural way the space and internal symmetries. We de- 
scribe the main features of the gauge principle relevant to the discrete 
and finite background. Assuming that continuous phenomena are ap- 
proximations of more fundamental discrete processes, we discuss ~ with 
the help of a simple illustration - relations between such processes and 
their continuous approximations. We propose an approach to introduce 
quantum structures in discrete systems, based on finite gauge groups. In 
this approach quantization can be interpreted as introduction of gauge 
connection of a special kind. We illustrate our approach to quantization 
by a simple model and suggest generalization of this model. One of the 
main tools for our study is a program written in C. 



1 Introduction 

In 1918 Hermann Wcyl - guided by the concept that the scale of length is 
arbitrary: if there is no fundamental length in Nature it does not matter what 
unit of length is used in measurements - conjectured that the scale can be taken 
in the form e'^'^^^ , i.e., it may vary from point to point in time and space. This idea 
failed in application to physics but gave rise to the concept of gauge invariance. 

Later - in 1929, after advent of Quantum Mechanics ~ Weyl (and also 
Vladimir Fock and Fritz London) replaced scale transformations e'^^^-' by ro- 
tations (phase transformations) e*'^'-^-' and derived electromagnetism from the 
gauge principle. 

In 1954 C.N. Yang and R. Mills extended the gauge principle to non-Abelian 
symmetries. Novif the gauge principle is recognized as one of the central principles 
in contemporary physics - in fact, all fundamental physical theories are gauge 
theories (for historical review see [T]). 

The lattice gauge theory was introduced by K.G. Wilson in 1974 as a practi- 
cal approach to the problems of strong interactions for which the standard per- 
turbative methods are inapplicable. This technique - based on approximation of 
space, or space-time, by some (usually hypercubic) lattice - was considered as an 



auxiliary computational method rather than a fundamental construction. The 
later mathematical generalizations established relations between lattice gauge 
theories and such topics as topological quantum field theory (TQFT), invariants 
of 3- and ^-manifolds, monoidal categories, Hopf algebras and quantum groups, 
quantum gravity etc., [2]. 

In view of their origin and applications, the above mentioned lattice gauge 
theories are not entirely discrete constructions. They involve continuous ingre- 
dients: gauge groups arc Lie groups, Lagrangians and observables are real or 
complex functions. Furthermore, the gauge groups of these theories are groups 
of internal symmetries and do not involve the lattice symmetries. It seems de- 
sirable to include the space symmetries into construction of gauge group, since: 
(a) the quantum statistics of particles is characterized by the rules describing 
their behavior under permutations of the points of space; (b) there exist gauge 
theories that deduce gravity by interpreting the space or space-time symmetries 
as gauge groups. 

In this paper we consider more radical version of discrete gauge invariance. 
All our manipulations including quantization remain within the framework of 
exact discrete mathematics requiring no more than the ring of algebraic integers 
(and sometimes the quotient field of this ring). Our study was carried out with 
the help of a program in C we are developing now. 

2 Discrete Dynamics 

Wc consider evolution in the discrete time ieZ = {...,— 1,0,1,...}. 

Let the space X be a finite set of points: X = . . . , xj^^}. This - primor- 
dially amorphous - set may possess some structure: some points may be "closer" 
to each other than others. A mathematical abstraction of such a structure is an 
abstract simplicial complex - a collection of subsets of X (simplices) such that 
any subset of a simplex is also simplex. Onc-dimensional complexes, i.e., graphs 
(or lattices), are sufficient to formulate a gauge theory. The symmetry group of 
the space X is the graph automorphism group G — Aut (X). 

Table [T] shows some lattices with their symmetries. We use these lattices 
in our computer experiments. In the table Nx and Ne are numbers of points 
and vertices in space X; the trivial, symmetric, cyclic, dihedral and alternating 
groups arc denoted by 1, Sym (n). Z„, D2„ and Alt(ri,), respectively; the signs x 
and XI denote direct and semidirect products, respectively. Note, that the lattice 
denoted as Toric square nxn'm the table has three times larger symmetry group 
at n = 4 than the general case formula predict^ 

Let each point x € X take values in some finite set of local states S = 
{(71, . . . , ctn^} possessing some symmetry group F < Sym (S) . Such groups are 
analogs of the "groups of internal symmetries" responsible for interactions in 



^ N. Vavilov pointed out to the author that tliis extra symmetry can be explained by 



Zg symmetry of the Dynkin diagram D4 ~ 




associated with the case n = 4. 
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physical gauge theories. The state of a system as a whole is a function a{x) G . 

Dynamics of the system is determined by some evolution rule connecting the 
current state of the system atix) with its prehistory at~i{x), at-2{x), ... . 
A typical form of evolution rule is evolution relation: 

R {<Jt {x) ■ at-i {x) , at-2 (a;) , . . .) C x 17^ x • • • . (1) 

Most commonly used in applications and convenient for study arc determinis- 
tic (or causal) dynamical systems. The current state of deterministic system is 
uniquely determined by its prehistory, i.e., relations like ([1]) are functional and 
can be written in the form 

at (x) = F {at-i {x) , at-2 {x) , . . .) . 

There arc two important special types of non- deterministic dynamical systems: 

— lattice models in statistical mechanics - special instances of Markov chains; 

— discrete quantum systems obtained from classical systems by identification 
of their states with basis elements of complex Hilbert spaces. 

For these systems transition from one state to any other is possible with some 
probability controlled by additional structures: real (for Markov chains) or com- 
plex (for quantum systems) weights assigned to state transitions. In this paper 
we restrict our attention to the case of discrete quantum systems. 

3 Unification of Space and Internal Symmetries 

Having the groups G and F acting on X and U, respectively, we can combine 
them into a single group W which acts on the states E-^ of the whole system. 
The group W can be identified, as a set, with the Cartesian product F^ (g) G, 
where F^ is the set of P-valued functions on X. That is, every element u G W 
can be represented in the form u — {a{x), a) , where a(x ) e F^ anda€ G. 

In physics, it is usually assumed that the space and internal symmetries are 
independent, i.e., W is the direct product F^ x G with actioro on S-^ and 
multiplication rule: 

a{x){a{x), a)^a-{x)a{x) action , 

{a{x), a)*{p{x), b) = {a (x) P (x) , ab) multiplication . (2) 

Another standard construction is the wreath product F Ix G having a structure 
of the semidirect product F-^ x G with action and multiplication 

(T{x){a{x), a) = a (^xa~^^ a (^xa~^) , 
{a{x), a)'f'{(3{x), b) = {a{x)(3{xa), ab) . (3) 

^ We write group actions on the right. This, more intuitive, convention is adopted 
in both GAP and MAGMA - the most widespread computer algebra systems with 
advanced facilities for computational group theory. 
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These examples are generalized by the following Statement: 
There are equivalence classes of split group extensions 1 W ^ G ^ 1 

determined by antihomomorphisms ^ : G G. The equivalence is described by 
arbitrary function k : G ^ G. The explicit formulas for main group operations 
— action on , multiplication and inversion — are 

a{x) (a (a;), a) 
(a(x), a) * iPix), h) 

{a{x), 

This statement follows from the general description of the structure of split 
extensions of a group Gi by a group Gq: all such extensions are determined by 
the homomorphisms from Gi to Aut (Go) (see, e.g., [3], P- 18). Specializing this 
description to the case when Go is the set of F-valued function on X and Gi 
acts on arguments of these functions we obtain our statement. The equivalence 
of extensions with the same antihomomorfism but with different functions k 
is expressed by the commutative diagram 

1 > > W > G > 1 



K 









1 > > w' > G > 1 

where the mapping K takes the form K : {a{x), a) i— > (a {xK(a)), a) . 

Note that the standard direct ^ and wreath ^ products are obtained from 
this general construction by choosing (/i(a) = l,K(a) = 1) and (/^(a) = a~^, 
K{a) = a~^), respectively. 

In our C program the group W is specified by two groups G and F and two 
functions /i(a) and K{a) implemented as arrays. It is convenient in computations 
to use the following specialization: /i(a) = a~™ and K{a) = . For such a choice 
formulas (l4][6]) take the form 

a{x) (a (x), a) 
{a{x), a) * (/5(x), b) 
{a{x), a) 

Here k is arbitrary integer, m = [direct product) or m = 1 [wreath product). 

4 Discrete Gauge Principle 

In fact, the gauge principle expresses the very general idea that any observable 
data can be presented in different "frames" at different points of space and 



= a [xiJL[a)) a[xK[a)) , (4) 
= [a{xK[ab)~^n[b)K[a)) (3{xK[aby\[b)), ab) , (5) 

= (a(^xK[a^^^ ^ fi[a)~^K[a)^ , '^^^^ • i^) 



= a [xa-"') a [xa'') , (8) 
= {a{x[ab)-''-"'a''+"') (3{x[ab)-''b''), ab) , (9) 
= (a[xa^''+"y\ a-i) . (10) 
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time, and there should be some way to compare these data. At the set-theoretic 
level, i.e., in the form suitable for both discrete and continuous cases, the main 
concepts of the gauge principle can be reduced to the following elements 

— a set X, space or space-time; 

— a set Z", local states; 

— the set S-^ of Z'-valued functions on X, the set of states of dynamical system; 

— a group W < Sym {S'^) acting on S'^ , symmetries of the system; 

— identification of data describing dynamical system with states from ZJ-^ 
makes sense only modulo symmetries from W; 

— having no a priori connection between data from S-^ at different points x 
and y in time and space we impose this connection (or parallel transport) 
explicitly as VF- valued functions on edges of abstract graph: 

Pix,y) G ^{y) = a{x)P{x,y) ; 

connection P{x,y) has obvious property P{y,x) = P{x,y)~^; 

— connection P{x, y) is called trivial if it can be expressed in terms of a function 
on vertices of the graph: P{x,y) = p{x)~^p{y), p{x),p{y) € W; 

— invariance with respect to gauge symmetries depending on time or space 
u{x),u{y) € W leads to transformation rule for connection 

Pix,y)^u{x)-'P{x,y)u{y); (11) 

— the curvature of connection P{x, y) is defined as the conjugacy class of the 
holonomy along a cycle of a graph: 

P{xi,X2, ...,Xk) = P{xi,X2)Pix2,X3) ■ ■ ■P{xk,Xi) 

(the conjugacy means P'{xi, . . . , Xk) ~ u^^P{xi, . . . , Xk)u for any u e W); 
the curvature of trivial connection is obviously trivial: P{xi, . . . ,Xk) = 1; 

— the gauge principle does not tell us anything about the evolution of the con- 
nection itself, so gauge invariant relation describing dynamics of connection 
{gauge field) should be added. 

Let us give two illustrations of how these concepts work in continuous case. 

Electrodynamics. Abelian prototype of all gauge theories. Here the 
set X is 4-dimensional Minkowski space with points x = (a;^) and the set of 
states is Hilbert space of complex scalar (Schrodinger equation) or spinor (Dirac 
equation) fields ipi^)- The symmetry group of the Lagrangians and physical 
observables is W = U(l). The elements of U(l) can be represented as e~*". 

Let us make these elements dependent on space-time and consider the parallel 
transport for two closely situated space-time points: 

P{x,x + Ax) = e-"'(^'^+'^^) . 



6 



Specializing transformation rule (jlip to this particular case 

P'{x, x + Ax) = e'"(^)p(a;, x + Z\x)e-*"(^+^^) , 
substituting approximations 

P{x, x + Ax)^ ^-ip(x.x+Ax) _ _ iA{x)Ax , 

P'{x, X + Ax) = e-'''(^'^+^^) « 1 - , 
e-'-^i'=+'^'=) « e-*"(^) (1 - i\/a{x)Ax) , 
and taking into account commutativity oiW ~ U(l) we obtain 

(J(X\X) 

A'{x) = Alx) + Vaix) or, in components, A'Jx) = AJx) + ^ . (12) 

The 1-form ^ taking values in the Lie algebra of U(l) and its differential F = 
(-P/iiv) = dA are identified with the electromagnetic vector potential and field 
strength, respectively. To provide the gauge invariance of the equations for field 
ip^x) we should replace partial by covariant derivatives 

df,^ D^^d^~ iA^{x) 

in those equations. 

Finally, evolution equations for the gauge field A{x) should be added. In the 
case of electromagnetics these are Maxwell's equations: 

d_F = first pair (13) 
d*-F = second pair. (14) 

Here * is the Hodge conjugation (Hodge star operator). Note that equation (jl4p 
corresponds to vacuum Maxwell's equations. In the presence of the current J 
the second pair takes the form * d * _F = J. Note also that the first pair is 
essentially a priori statement, it reflects simply the fact that F, by definition, is 
the differential of an exterior form. 



Non-Abelian gauge theories in continuous space-time. Only minor mod- 
ifications are needed for the case of non-Abelian Lie group W. Again expansion 
of the VF- valued parallel transport for two close space-time points x and x + Ax 
with taking into account that P{x,x) = 1 leads to introducing of a Lie algebra 
valued 1-form A = (Af^) : 

P{X,X + Ax) + Af,{x)Ax^' . 

Infinitesimal manipulations with formula pip 

u{x)~'^P{x,x+Ax)u{x+Ax) u{x)-^{l + Af,{x)Ax'') (u{x) + ^^^AxA 
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lead to the following transformation rule 



A'^{^)^u{x)-'A,{x)u{x)^u{x)-'^ . (15) 
The curvature 2-forni 

F ^dA+[AhA] 

is interpreted as physical strength field. In particular, the trivial connection 

, 9mo(x) 



Ai_,{x) = ito(x)" 



^x^' 



is flat, i.e., its curvature F = 0. 

There are different approaches to construct dynamical equations for gauge 
fields [2]. The most important example is Yang-Mills theory based on the La- 
grangian 

Lym = Tr [F a Tk-F] . 
The Yang-Mills equations of motion read 

dF + A F] = , (16) 
d*F+ [A A*F] = . (17) 

Here again equation (jl6p is a priori statement called Bianci identity. Note that 
Maxwell's equations are a special case of Yang-Mills equations. 

It is instructive to see what the Yang-Mills Lagrangian looks like in the dis- 
crete approximation. Replacing the Minkowski space X by a hypercubic lattice 
one can see that the discrete version of Lym is proportional to cr (7/), where 
the summation is over all faces of a hypercubic constituent of the lattice; 

(7 = 2 dim U - {xu + Xc/t ) ; 

Xu and Xc/t are characters of the fundamental representation U of the gauge 
group and its dual representation, respectively; 7/ is the gauge group holonomy 
around the face /. 

The Yang-Mills theory uses Hodge operation converting fc-forms to {n — k)- 
forms in ri-dimensional space with metric g^^,^. In topological applications so- 
called BF theory plays an important role since it does not involve a metric. 
In this theory, an additional dynamical field B is introduced. The Lie algebra 
valued {n — 2)-form B and the 2-form F are combined into the Lagrangian 
LBF = Tr[BAF] . 



5 Quantization Based on Finite Group 

Quantization is a procedure for recovering a more fundamental quantum theory 
from its classical approximation. Both Lagrangian and Hamiltonian formulations 
of classical mechanics are based on the principle of least action which looks a 
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bit mysterious: a particle moving from one point to another "knows" in advance 
where it is going to arrive. Feynman's path integral quantization [3] eliminates 
this apparent teleology in a quite natural way: classical trajectories correspond 
to the dominating (in the path integral) part of all possible trajectories. 

Of course, recovering a theory from its approximation can not be performed 
uniquely. Moreover, discrepancies between a theory and its approximation may 
be essential. To illustrate, let us compare a simple discrete process with its 
approximation by continuous physical law. 



5.1 Heat Equation from Bernoulli Trials 

Let us consider a sequence of Bernoulli trials. The probability of a separate 
sequence is described by the binomial distribution 

P(n_,n+ = 1 j— p_ . (18) 

Here {— , +} are possible outcomes of a single trial; p_, p+ are probabilities 
(p_ +P+ = 1) and n_, n_|_ are numbers of the outcomes. 

Applying Stirling's approximation to (jlSp and introducing new variables 
X ~ 71+ — n_, t = 71- + 7i+, V = p- — p+ — let us call them "space", "time" 
and "velocity" , respectively — we obtain 

P{x,t)^P{x,t)= W£-exp <;---( ^== ) 1 . (19) 




This is the fundamental solution of the heat (also known as diffusion or Fokker- 
Planck) equation: 

dP{x,t) dP{x,t) _ {l^v^)d^P{x,t) 

~^r~ + ''^^ - ~^ ■ 

Note that expression (fig]) contains "relativistic" fragment ^ ~ ^'^ due to 

Vl — f 

the velocity limits —1 < w < 1 in our model. Note also that at \v\ = 1 equation 
(|20p reduces to the wave equation 

dP{x,t) , dP{x,t) _^ 
dt ^ dx " " ■ ' 

Now let us set a problem as is typical in mechanics: find extremal trajectories 
connecting two fixed points (0,0) and (X, T). We adopt here the search of tra- 
jectories with maximum probability as a version of the "least action principle". 
The probability of trajectory passing through some intermediate point {x, t) is 
the following conditional probability 

_P{x,t)P{X~x,T~t) 

^(0,0)^(x,t)-,(X,T) — 



P{X, T) 

t\{T-t)\{^^\UI±2L\\ 



2 )■ \ 2 



g-a; '^! (*+£)] ^ T-t _ X-x y^ ^ T-t _|_ X-x yrp\ 



(22) 
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The conditional probability computed for approximation (|19[) takes the form 



~ _ ^ / (Xt-xTf 

P(o,o)^(.,*)^(^,T) - -j======== exp j-2(i_„2)ir(r_i) 

One can see essential differences between and 

— exact probabilities ()22p do noi depend on the velocity w (or on the probabil- 
ities p+ of a single trial), whereas (j23p contains artificial dependence, 

— it is easy to check that expression p2l) allows many trajectories with the 
same maximum probability, whereas extremals of (|23p are deterministic tra- 
jectories, namely, straight lines x = -^i. 

These artifacts show that an important guiding principle of quantization — 
correspondence with classical limit — may not be quite reliable. 

5.2 Gauge Connection and Quantization 

The Aharonov-Bohm effect (Fig. [T|) is one of the most striking illustrations of 
interplay between quantum behavior and gauge connection. Charged particles 
moving through the region containing perfectly shielded thin solenoid produce 
different interference patterns on a screen depending on whether the solenoid is 
turned on or off. There is no electromagnetic force acting on the particles, but 
working solenoid produces U(l)-connection adding or subtracting phases of the 
particles and, thus, changing the interference pattern. 




Fig. 1. Aharonov-Bohm effect. Magnetic flux is confined within the perfectly 
shielded solenoid; interference pattern is shifted in spite of absence of electro- 
magnetic forces acting on the particles. 



In the discrete time Feynman's path amplitude decomposes into the product 
of elements of the group U(l) (or, more precisely, elements of the fundamental 
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representation of U(l)): 

exp(,-S,.cxp(./«.)^c-o,.^„„-.-»^, ^ ,24) 

By the notation it-i.t we emphasize that the Lagrangian is in fact a function 
defined on pairs of points (graph edges) — this is compatible with physics where 
the typical Lagrangians are determined by the first order derivatives. Thus, the 
expression P{t — l,t) = g^Lt-i^t ^ U(l) can be interpreted as U(l)-parallel 
transport. A natural generalization of this is to suppose that: 

— the group U(l) can be replaced by some other group F, 

— unitary representation p{r) may have dimensionality different from 1. 

We can introduce quantum mechanical description of a discrete system in- 
terpreting states a & E as basis elements of a Hilbcrt space This allows to 
describe statistics of observations of a's in terms of the inner product in ^ . 

Now let us replace expression for Feynman's path amplitude by the 
following parallel transport along the path 

^p(r) = P (ao,i) • ■ ■ P • ■ • P ("t-i^t) • 

Here at-i.t are elements of a finite group F - we shall call F quantizing group - 
and p is an unitary representation of F on the space ^ . 

Let us recall main properties of linear representations of finite groups [5] . 

— First of all, any linear representation of finite group is equivalent to unitary. 

— Any unitary representation p is determined uniquely (up to isomorphism) 
by its character defined as Xp{ot) = Tr/5(a), a € F. 

— All values of Xp ^-nd eigenvalues of p are elements of the ring A of algebraic 
integers, moreover the eigenvalues arc roots of unity. Recall that the ring A 
consists of the roots of monic polynomials with integer coefficients j^. 

— If all different irreducible representations of F are pi,--- ,pi,--- , ph and 
di = dim Pi, M = \F\ then 

h 

^df = M and any divides M : d^\ M. 

i-l 

— Any function (p{a) depending only on conjugacy classes of F, i.e., 

ip (^f3~^al3^ = ip (a), is linear combination of characters Xpi ? ■ ' ■ > Xph ■ 
Such functions are called central or class functions. 

If the group F consists of M elements 70, . . . ,7a/-i and is the number of 
paths with the "phase" ^ = pijk) at the point of observation {x,t), then the 

amplitude at this point is A = ^ rikp{'yk) V': where ip The square of the 

fc=0 

amplitude (i.e., probability after appropriate normalization) can be written as 

A/-1 

{A^\Ai,) = + J2 ^^"fc(^|p(7rV^) + pU7rSfc)|v^) , (25) 

i<k 
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or, after collecting like terms, as 

M-l 

(AV-lvl^) ^ iVfc(no,--- ,nM-i)(V'|/3(7fe)+P^(7fc)|V') , (26) 

A;=0 

where (no, ■ • • , um-i) are quadratic polynomials with integer coefficients and 
arguments. Thus, algebraic integers are sufficient for all our computations except 
for normalization of probabilities requiring the quotient field of the ring A. 

5.3 Simple Model Inspired by Free Particle 

In quantum mechanics - as is clear from the never vanishing expression exp (-^iS*) 
for the path amplitude - transitions from one to any other state are possible in 
principle. However, we shall consider computationally more tractable models 
with restricted sets of possible transitions. 

Let us consider quantization of a free particle moving in one dimension. Such 

a particle is described by the Lagrangian L = ■ Keeping only transitions to 
the closest points in the discrctized space we come to the following rule for the 
one-time-step transition amplitudes 

i_ m{(x+l)-x}'^ ■ m 

" ' e 



1 

i— 
e 2h. 

That is, we have evolution rule as an U(l)-valued function R defined on pairs of 
points (graph edges). Symbolically: 

i?(a; ^ .t) = 1 e U(l) , 

■ m 

R{x^x-l)^R{x->x + l)^w^e^2h e U(l). (27) 

Now let us assume that w in (j27p is an clement of some representation of a 
finite group: w ~ p (a) , a G P = {70 = 1, . . . , 7m-i}- Rearranging multinomial 
coefficients — trinomial in this concrete case — it is not difficult to write the 
sum amplitude over all paths of the form (0, 0) — > {x,t) 

Note that x must lie in the limits determined hy t: x € [—t, t] . 

One of the most expressive peculiarities of quantum-mechanical behavior is 
the destructive interference — cancellation of non-zero amplitudes attached to 
different paths converging to the same point. By construction, the sum of ampli- 
tudes in our model is a function A(w) depending on distribution of sources of the 
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(28) 



particles, their initial phases, gauge fields acting along the paths, restrictions - 
like, e.g., "slits" ~ imposed on possible paths, etc. In the case of one-dimensional 
representation the function A^w) is a polynomial with algebraic integer coeffi- 
cients and w is a root of unity. Thus, the condition for destructive interference can 
be expressed by the system of polynomial equations: A{w) = and iv'^^ = 1. For 
concrctcncss let us consider the cyclic group F = Zm = {70, ■ • • , 7^, • ■ ■ , 7m-i}- 
Any of its M irreducible representations takes the form pijk) — , where w 
is one of the A/th roots of unity. For simplicity let w be the primitive root: 
w = e^'^*/^'^. Fig. [2] shows all possible transitions from the point x in three time 
steps with their amplitudes. 




x±2 



Number of paths 7 6 3 

Amplitude Af, = I + Gw'^ AL = 3w + 3w'^ At, = 3uP 



Fig. 2. Amplitudes for all possible paths in three time steps. 



We see that the polynomial Aj.-^ = 3w + 3w'^ — 3w {w^ + l) contains the 
cyclotomic polynomial (^4(111) = + 1 as a factor. The smallest group associated 
to ^^{w) — and hence providing the destructive interference — is Z4. Thus, Z4 
is the natural quantizing group for the model under consideration. 

Fig. [3] shows interference patterns — normalized squared amplitudes ("prob- 
abilities") — from two sources placed in the positions x ~ —4 and x = A for 20 
time steps. The upper and lower graphs show interference pattern when sources 
are in the same {Acp = 0) and in the opposite (Acj) = tt) phases, respectively. 



5.4 Generalization: Local Quantum Model on Regular Graph 

The above model — with quantum transitions allowed only within the neigh- 
borhood of a vertex of a 1-dimensional lattice — can easily be generalized to 
arbitrary regular graph. Our definition of local quantum model on k-valent graph 
uncludes the following: 

1. Space X = {xi, ■ ■ ■ ,xn} is a fc— valent graph. 

2. Set of local transitions Ei = {eo,i, ei^i, • • • , ek,i} is the set of k adjacent to the 
vertex Xi edges em,i = (xi Xm,i) completed by the edge eo,i = {xi Xi). 

3. We assume that the space symmetry group G = Aut {X) acts transitively 
on the set {Ei, • • • , E^}. 

4. Gi = Stabc {xi) < G is the stabilizer of Xi {g € Gi means xtg = .t,). 

5. f2i = {a;o,i, • • • ,uJh,i} is the set of orbits of Gi on Ei. 
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Fig. 3. Group Z4. Interference from two sources. Number of time steps T = 20. 
Source positions are -4 and 4. Phase differences A(f) = 04 — between sources 
are and tt. 



6. Quantizing group _r is a finite group: F = {70, • ■ ■ ,7a/-i}. 

7. Evolution rule i? is a function on Ei with vahies in some representation p (F). 
The rule R prescribes p (P)-wcights to the one-time-step transitions from Xi 
to elements of the neighborhood of .t,; . From the symmetry considerations R 
must be a function on orbits from i?.;, i.e., R {em,ig) ~ R {em,i) for g G G;. 

To illustrate these constructions, let us consider the local quantum model on the 
graph of buckyball. The incarnations of this 3-valent graph include in particular: 

- the Caley graph of the icosahedral group Alt (5) (in mathematics); 

- the molecule Ceo (in carbon chemistry) . 



Here the space X = {xi, • • • ,xqq} has the shape 




and its symmetry 



group is G = Aut {X) = Z2 x Alt (5). The set of local transitions takes the form 
Ei = {eoj, eij, e2j, 63,,;}, where eo,i = {x, Xi), eij = {xi xi^i), 



62/ 



^2,i), 63. i = {xi x^j) in accordance with 




The stabilizer of Xi is Gi ~ Stabc (xi) ~ Z2. The set of orbits of Gi on Ei con- 
tains 3 orbits: f2i = {wq,, = {eo,i},wi,, = {ei,i, 62,^} , a;2,i = {e3,i}}, i.e., the sta- 
bilizer does not move the edges {xi Xi) and {xi —>■ x^^i) and swaps {xi Xi^i) 
and {xi — > X2.i) ■ This asymmetry results from different roles the edges play in 
the structure of the buckyball: {xi xi.i) and {xi — > X2,i) are edges of a pen- 
tagon adjacent to Xi , whereas {xi X3^i) separates two hexagons; in the carbon 
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molecule Cgo the edge {xi — *■ x^^i) corresponds to the double bond, whereas oth- 
ers are the single bonds. 

The evolution rule takes the form: 

R {xi Xi) = p (ao) , 

R {xi xi^i) = R{xi X2,i) = p (ai) , 

R{xi X3^i) = p{a2) , 

where ao, ai, Q!2 G r . If we take a one-dimensional representation and move ao 
- using gauge invariance - to the identity element of F , we sec that the rule R 
depends on v ~ p{oLi) and w = p{a2). Thus, the amplitudes in the quantum 
model on the buckyball take the form A{v, w) depending on two roots of unity. 

6 Conclusion 

Extraordinary success of gauge theories in fundamental physics suggests that 
the gauge principle may be useful in theory and applications of discrete dynam- 
ical systems also. Furthermore, discrete and finite background allowing compre- 
hensive study ~ especially with the help of computer algebra and methods of 
computational group theory - may lead to deeper understanding of the gauge 
principle itself and its connection with the quantum behavior. To study more 
complicated models we are developing the C program. 

Acknowledgments 

The author thanks Laurent Bartholdi, Vladimir Gerdt and Nikolai Vavilov for 
useful remarks and comments. This work was supported in part by the grants 
07-01-00660 from the Russian Foundation for Basic Research and 1027.2008.2 
from the Ministry of Education and Science of the Russian Federation. 

References 

1. O'Raifeartaigh L., Straumann N. Gauge theory: Historical Origins and Some Mod- 
ern Developments. Reviews of Modem Physics, 72, No. 1 (2000) 1-23 

2. Oeckl R. Discrete Gauge Theory (From Lattices to TQPT). Imperial College Press, 
London (2005) 

3. Kirillov A. A. Elements of the Theory of Representations. Springer- Verlag, Berlin- 
New York (1976) 

4. Feynman R.P., Hibbs A.R. Quantum Mechanics and Path Integrals. McGraw-Hill 
(1965) 

5. Serre J. -P. Linear Representations of Finite Groups. Springer- Verlag (1977) 



15 



